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ABSTRACT 

We follow the changes in the structure of the accretion disk of the dwarf nova V2051 Oph 
along two separate outbursts in order to investigate the causes of its recurrent outbursts. We 
apply eclipse mapping techniques to a set of light curves covering a normal (July 2000) and a 
low-amplitude (August 2002) outburst to derive maps of the disk surface brightness distribution 
at different phases along the outburst cycles. The sequence of eclipse maps of the 2000 July 
outburst reveal that the disk shrinks at outburst onset while an uneclipsed component of 13 per 
cent of the total light develops. The derived radial intensity distributions suggest the presence 
of an outward-moving heating wave during rise and of an inward-moving cooling wave during 
decline. The inferred speed of the outward-moving heating wave is ~ 1.6 kms~^, while the speed 
of the cooling wave is a fraction of that. A comparison of the measured cooling wave velocity on 
consecutive nights indicates that the cooling wave accelerates as it travels towards disk center, 
in contradiction with the prediction of the disk instability model. From the inferred speed of the 
heating wave we derive a viscosity parameter ahot — 0.13, comparable to the measured viscosity 
parameter in quiescence. The 2002 August outburst had lower amplitude {AB 0.8 mag) and 
the disk at outburst maximum was smaller than on 2000 July. For an assumed distance of 92 pc, 
we find that along both outbursts the disk brightness temperatures remain below the minimum 
expected according to the disk instability model. The results suggest that the outbursts of V2051 
Oph are caused by bursts of increased mass transfer from the mass-donor star. 

Subject headings: stars: dwarf novae, cataclysmic variables — stars: individual (V2051 Ophiuchi) — 
accretion, accretion disks 

1. Introduction 
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to a white dwarf companion via an accretion disk 
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or column (Warner 1995). The subclass of dwarf 
novae comprises low-mass transfer CVs in which 

mass is fed to a weakly magnetic ( B < 10®G) 

white dwarf. These binaries show recurrent out- 
bursts of up to 5 magnitudes on timcscales of 
weeks to months as a consequence of a sudden in- 
crease in mass inflow in the accretion disk. 

There are two competing models to explain the 
cause of the sudden increase in mass accretion. In 
the disk- instability model (DIM), matter is trans- 
ferred at essentially constant rate to a low vis- 
cosity disk and accumulates in an annulus until a 
thermal-viscous instability switches the disk to a 
high viscosity regime and the gas diffuses rapidly 
inward and onto the white dwarf. In this case, the 
outburst starts at the radial position where the 
thermal instability first occurs and propagates as 
a heating wave on the rise and as a cooling wave 
along the decline back to quiescence (Hameury 
et al. 1998; Lasota 2001). In the mass transfer in- 
stability model (MTIM), the outburst is the time- 
dependent response of a high viscosity accretion 
disk to a burst of enhanced mass transferred from 
the secondary star (Bath 1975). This model pre- 
dicts that the disk shrinks at the onset of the out- 
burst in response to the sudden addition of matter 
with low angular momentum, and that the disk 
viscosity in quiescence and in outburst are simi- 
lar. DIM predicts no reduction of disk radius at 
outburst start but demands that the viscosity pa- 
rameter in quiescence be 5-10 times smaller than 
the viscosity in outburst. The thermal limit cycle 
of DIM also predicts that the outbursting parts 
of the accretion disk should be hotter than a crit- 
ical temperature rcfi(crit) (Warner 1995). Fur- 
thermore, the inward cooling wave is expected to 
decelerate as it travels towards disk center (Menou 
et al. 1999). Thus, there arc several predictions 
that can be tested from observations of accretion 
disks through outburst cycles with the aid of in- 
direct imaging techniques such as eclipse mapping 
(Home 1985). 

The interest in testing both models against ob- 
scirvations has somehow reduced over the last two 
decades as a consequence of a wide acceptance of 
DIM as the correct explanation. Two arguments 
have been key in setting the dominance of DIM. 
They are based on the assumption that the matter 
transferred from the secondary star is deposited at 



the disk rim, where the gas stream hits the accre- 
tion disk to form a bright spot (BS). Given this 
assumption, one may predict that (i) a burst of 
enhanced mass transfer rate would inevitably lead 
to an increase in the luminosity of the BS and that 
(ii) an MTIM-driven outburst can only lead to 
outside-in outbursts because the additional matter 
is always deposited at the disk rim. The existence 
of inside-out outbursts and the lack of observa- 
tional support for an increase in BS luminosity at 
outburst onset seem to argue against MTIM. How- 
ever, one should note that both arguments (and 
the reasoning to drop MTIM) fall apart if the as- 
sumption upon which they rely is incorrect. As 
we shall see later (Sect. 5), this may indeed be the 
case. 

V2051 Oph is an ultra-short-period CV {Porb < 
2 hr), almost all of which are either polars with 
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highly magnetic white dwarfs 

members of the SU UMa subclass of superout- 
bursting dwarf novae (Warner 1995). The binary 
was discovered by Sanduleak (1972) and since then 
several classifications were proposed. Warner & 
O'Donoghue (1987) listed some characteristics of 
the object that distinguish it among other CVs: 
with its Porb = 90 min it has one of the shortest 
periods known; it is one of the few short-period 
CVs to show deep eclipses {B ~ 2.5 mag); unlike 
other eclipsing systems, it possesses strong flicker- 
ing (random brightness fluctuations of 0.1-1 mag- 
nitudes on timescales from seconds to minutes) 
and flare activity which producx; a wide assortment 
of eclipse morphologies (Warner & Cropper 1983). 
V2051 Oph was proposed to be a low-field polar 
system by Warner & O'Donoghue (1987) based 
on the interpretation of their eclipse maps and on 
the observation of a 42-s oscillations in the optical 
during outburst, reminiscent of the rapid oscilla- 
tion seen in polars. A consensus was reached with 
the observation of a superoutburst during which 
superhumps were detected (Kiyota & Kato 1998, 
Vrielmann & Offutt 2003), establishing its classi- 
fication as an SU UMa-type dwarf novae. 

As part of a long-time campaign to study the 
flickering properties of a sample of eclipsing CVs 
(e.g., Baptista & Bortoletto 2004), we happened 
to observe V2051 Oph in outburst on two occa- 
sions. This paper reports the analysis of the light 
curves of V2051 Oph collected along these two out- 
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bursts with eclipse mapping techniques. The re- 
sulting maps deliver the disk brightness distribu- 
tion in quiescence, on the rise to maximum, dur- 
ing maximum light and in the decline phase. Sec- 
tion 2 describes the observations and Sect. 3 gives 
the details of the data analysis. The results are 
presented in Sect. 4 and discussed in Sect. 5 in a 
comparison with predictions from the dwarf nova 
outburst models. 

2. The observations 

Time-series of high-speed CCD photometry of 
V2051 Oph were obtained with a CCD camera 
(385 X 578 pixels, 0.58 arcsec/pixel) attached to 
the 1.6m telescope of Laboratorio Nacional de As- 
troffsica (LNA), in Itajuba, Brazil, between 2000 
July 28 and August 2 in B and V bands, and on 
2002 August 4, 5 and 7 in the B band. The ob- 
servations are summarized in Table [T] The sec- 
ond column lists the starting time of the observa- 
tions in Julian Date and the third column shows 
the time resolution in seconds (At). The fifth col- 
umn lists the eclipse cycle number (E); the num- 
bers in parenthesis indicate observations that, be- 
cause of interruptions caused by, p. ex., clouds, do 
not cover the eclipse itself. The sixth column in- 
dicates the binary phases covered by the obser- 
vations (running from —0.5 to -1-0.5 with eclipse 
center at phase 0.0), the seventh column lists the 
number of points in the light curve (Np). The 
eighth column gives an estimate of the quality of 
the observations. The seeing ranged from 1.0 to 
2.2 arcsec. Part of the data (quiescence) was pre- 
sented in Baptista & Bortoletto (2004), but is re- 
peated here for completeness. All light curves were 
obtained with the same instrumental set and tele- 
scope, which ensures a high degree of uniformity 
to the data set. 

Data reduction procedures included bias sub- 
traction, flat-field correction, cosmic rays removal 
and aperture photometry extraction. Time series 
were constructed by computing the magnitude dif- 
ference between the variable and a reference com- 
parison star with scripts based on the aperture 
photometry routines of the APPHOT/IRAF pack- 
ag^. Light curves of other comparison stars in 
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the field were also computed in order to check the 
quality of the night and the internal consistency 
and stability of the photometry over the time span 
of the observations. 

We estimate that the absolute photometric ac- 
curacy of these observations is about 10 per cent. 
On the other hand, the analysis of the relative flux 
of the comparison star of all observations indicates 
that the internal error of the photometry is less 
than 2 per cent. The individual light curves have 
typical signal-to-noise ratios oi S/N = 40 — 50 out- 
of-eclipse and S/N = 10 — 20 at mid-eclipse. Ad- 
ditional details of the data reduction procedures 
are given in Baptista & Bortoletto (2004). 

The historical hght curve of V2051 Oph E 
around the epoch of both outbursts is shown in 
Fig. HI The star went in outburst between 2000 
July 30 and 31 (about JD 2451758) [Fig.IB top 
panel]. This was a normal outburst, with an am- 
plitude of 2 magnitudes and a time span of 2 days. 
Our observations cover three days before maxi- 
mum, the outburst maximum, and two days along 
the decline phase. The star was also seen in out- 
burst on 2002 July 9-10 (Fig.[Tl bottom panel). We 
did not observe this outburst. Our observations 
at this epoch cover the following, low-amplitude 
(~ 0.8 mag) and short outburst on 2002 August 
5 (JD 2452493). The runs on 2002 Aug 4 and 7 
frame, respectively, the day before outburst and a 
late decline stage when the star was almost back 
to its previous quiescent brightness level. A simi- 
lar, low-amplitude outburst occurred 20 days after 
the one recorded in our data. 

Average out-of-eclipse B and V magnitudes 
measured from our light curves are superimposed 
on the AAVSO historical data (Fig.[Il). The B- 
band magnitudes were displaced vertically by —0.5 
mag for visualization purposes. Only B-band data 
is available for the 2002 outburst. The V-band 
magnitudes are consistent with the visual obser- 
vations. However, the B-band magnitudes are 
systematically fainter than both the V-band and 
visual contemporary measurements, implying a 
color index {B — V) ~ -1-0.5 mag. This is no- 
ticeably redder than the color of V2051 Oph previ- 
ously reported in the literature, —0.1 < {B — V) < 
+0.2 (Bruch 1983; Vogt 1983), and may be asso- 
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ciated to the extended state of lower brightness 
the object has gone through in 1996 (see Bap- 
tista et al. 1998a). The color index changed from 
(B-V)^ +0.4 mag to (B-V) ^ +0.6 mag from 
quiescence to outburst maximum in 2000. 
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Fig. 1. — Historical Hght curve of V2051 Oph 
around the epoch of the two outbursts, con- 
structed from observations made by the AAVSO 
(crosses). Average out-of-eclipse B and V mag- 
nitudes from our data are indicated respectively 
by filled squares and open symbols. The B band 
magnitudes were displaced vertically by —0.5 mag 
for visualization purposes. 



3. Data analysis 

3.1. Light curve construction 

The B- and V-band individual light curves of 
V2051 Oph are shown in Fig. [2l with arbitrary 
adjustments in the x-axis separation and V-band 
fluxes for visualization purposes. 

The data were grouped per passband and out- 
burst stage and combined to produce light curves 
of improved S/N and reduced flickering influ- 
ence. The B-band data from 2000 July 28 and 
29 were combined into a single average light curve 
representative of the quiescence state (hereafter 
Jul2829) because they are at the same brightness 
level and there is no perceptible difference in light 
curve morphology from one night to the other. 
The B-band data from 2000 July 30 (JulSOB) was 
treated separately as it shows a decrease in bright- 
ness with respect to the data from the previous 
two nights as well as changes in eclipse shape 
(Fig. [3]). There is a continuous decrease in bright- 
ness along the observations on 2000 Jul 31 and 
Aug 1. The individual light curves of those nights 
were scaled (by factors < 10 per cent) to a mean, 
common out-of-eclipse flux level before combining 
them. For both nights it was possible to obtain a 
good match of the eclipse shape and out-of-eclipse 
level for all combining curves. 

For each light curve, we divided the data into 
phase bins of 0.003 cycle and computed the me- 
dian flux at each bin to reduce the influence of 
flickering. The median of the absolute deviations 
with respect to the median flux was taken as the 
corresponding uncertainty for each bin. The light 
curves were phase-folded according to the linear 
ephemeris (Baptista et al. 2003), 

Trrnd = BJDD 2443 245.97752+0.062 427 8634x£;, 

a) 

where BJDD denotes Barycentric Dynamical 
Time. A small phase corrections of +0.0061 cycle 
was further applied to the 2002 data to remove the 
long-term cyclical period changes (Baptista et al. 
2003) and to make the center of the white dwarf 
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Fig. 2.— Light curves of V2051 Oph through the 2000 July (top) and 2002 August (bottom) outbursts. The 
separation of data from consecutive nights is arbitrarily reduced to 0.3 d for visualization purposes. Labels 
above the data of each night indicate the date. The V-band light curves are shown with gray symbols and 
are indicated by a small 'v' label above the data. The V-band fluxes are arbitrarily scaled by a factor 0.75 
also for visualization purposes. 
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eclipse coincident with phase zero. 

The white dwarf (WD) and bright spot (BS) 
echpses are seen as sharp changes in the slope of 
the eclipse shape in quiescence light curves (e.g., 
Fig. [3]) . We separated the contribution of the WD 
from the Jul2829, JulSOB and Aug04 average light 
curves with a light-curve decomposition technique 
(Wood et al. 1985) in order to compute eclipse 
maps of only the accretion disk. A zoom of the 
Jul2829 and JulSOB light curves around eclipse 
is shown in Fig. [3l The original light curve is 
smoothed with a median filter of width 0.006 cycle 
(the phase width of the WD ingress/egress fea- 
ture) and its numerical derivative is calculated. 
The derivative curve is smoothed with the same 
median filter as above to reduce noise and improve 
the detection of the WD and BS features. The WD 
and BS ingress/egress features are seen as those 
intervals where the derivative is significantly dif- 
ferent from zero. A spline function is fitted to 
the remaining regions in the derivative to remove 
the contribution of the slowly varying eclipse of 
the extended disk, and estimates of the WD flux 
are obtained by integrating the spline-subtracted 
derivative curve at ingress and egress. The light 
curve of the WD is reconstructed by assuming that 
its flux is zero between ingress and egress and 
constant outside eclipse. The WD contribution 
is then removed from the data by subtracting the 
reconstructed WD curve from the original light 
curve. The eclipse maps of the Jul2829, JulSOB 
and Aug04 data shown in Sect. 14.21 were computed 
using WD-subtracted light curves. 

A close inspection on the light curves in Fig. [3] 
reveal that, while the phase width of the WD 
ingress/egress feature is the same on both nights, 
the integrated WD flux is reduced by ~ 25 per 
cent on JulSOB. This could be caused by an in- 
creased optical depth of the inner disk gas or by a 
thicker disk rim on that night, resulting in a lower 
visible fraction of the WD surface. 

The average light curves are shown in the left 
panels of Figs. H [5] and El The Jul2829 and the 
Aug05 light curves have flares at ~ —0.11 cycle. 
The error bars around these fiares were artificially 
increased to minimize their influence on the eclipse 
mapping modeling. 



3.2. Eclipse mapping 

The eclipse-mapping method (Horne 1985) is 
an inversion technique that uses the information 
in the eclipse shape to reconstruct the disk sur- 
face brightness distribution. The reader is referred 
to Baptista (2001) for a review on the subject. 
Eclipse mapping techniques (Baptista & Steiner 
1993) were applied to the average light curves of 
Sect. 13.11 to solve for a map of the disk surface 
brightness plus the flux of an additional uneclipsed 
component in each case. The uneclipsed compo- 
nent accounts for all light that is not contained 
in the eclipse map (i.e., light from the secondary 
star, a vertically extended disk wind, or both). 
The reader is referred to Rutten et al. (1992a) 
and Baptista et al. (1996) for the details of and 
tests with the uneclipsed component. 

All variations in a light curve in the standard 
eclipse mapping method are interpreted as being 
caused by the changing occultation of the emitting 
region by the secondary star. Thus, any out-of- 
eclipse brightness change (e.g., the orbital hump 
caused by anisotropic emission from the BS) has 
to be removed before applying the technique to 
a light curve. This was done by fitting a spline 
function to the phases outside eclipse, dividing the 
light curve by the fitted spline, and scaling the re- 
sult to the spline function value at phase zero in 
each case. This procedures removes orbital mod- 
ulations with minimal effects on the eclipse shape 
itself. 

Our eclipse map is a flat Cartesian grid of 
75 X 75 pixels centered on the primary star with 
side 2Rli (where Rli is the distance from the disk 
center to the inner Lagrangian point LI). The 
eclipse geometry is defined by the mass ratio q 
and the inclination i, and the scale of the map 
is set by Rli- We adopted the values of Bap- 
tista et al. (1998a), Rli = 0.422 i?Q, q = 0.19 
and i — 83°, which correspond to a white dwarf 
eclipse width of A(p ~ 0.0662 cycle. This combi- 
nation of parameters ensures that the white dwarf 
is at the center of the map. The reconstruc- 
tions were performed with a polar gaussian de- 
faiflt function (Rutten et al. 1992a) with radial 
blur width Ar — 0.0266 Rli and azimuthal blur 
width A9 = 30°. The reconstructions reached a 
final near or equal 1 for all light curves. 

The statistical uncertainties in the eclipse maps 
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Fig. 3. — Left: The Jul2829 (top) and Jul30B (bottom) light curves and corresponding derivative curves at 
a phase resolution of 0.001 cycle. The median filtered version of the light curve is shown as a solid grey line 
in each case. Vertical dotted lines indicate the contact phases of the WD and vertical dashed lines mark the 
mid-ingress {(jibi) and mid-egress {<j>be) phases of the BS. Right: the binary geometry for q = 0.19. Dotted 
lines depict the primary Roche lobe and the ballistic stream trajectory. The WD is shown as a cross; the BS 
positions are marked by diamonds and the corresponding disk sizes are indicated by circles of radii 0.48 Rli 
(dashed line) and 0.38 Rli (solid line), respectively for Jul2829 and Jul30B. 
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are estimated with a Monte Carlo procedure (see, 
e.g., Rutten et al. 1992a). For a given input data 
curve, a set of 20 artificial light curves are gener- 
ated in which the data points are independently 
and randomly varied according to a Gaussian dis- 
tribution with standard deviation equal to the un- 
certainty at that point. The artificial curves are 
fitted with the eclipse-mapping algorithm to pro- 
duce a set of randomized eclipse maps. These are 
combined to produce an average map and a map 
of the residuals with respect to the average, which 
yields the statistical uncertainty at each pixel. A 
map of the statistical significance (or the inverse of 
the relative error) is obtained by dividing the true 
eclipse map by the map of the standard deviations. 
The uncertainties obtained with these procedures 
are also used to estimate the errors in the derived 
radial intensity and temperature distributions. 

4. Results 

4.1. Bright spot phases and disk radius 
changes 

We used the derivative technique of Sect. l3.1l to 
measure the mid- ingress/egress phases of the BS 
(06i,06e) and to estimate the disk radius under 
the assumption that the BS is located where the 
stream of transferred matter hits the edge of the 
accretion disk. and are taken as the phases 
of, respectively, minimum and maximum of the 
features corresponding to BS ingress/egress in the 
smoothed and spline-subtracted derivative of the 
light curve. 

BS mid- ingress/egress phases for the light 
curves Jul2829 and JulSOB are shown as ver- 
tical dashed lines in the left-hand panels of 
Fig. H We find (f)bt = -0.0170 ± 0.0015 and 
(t)be = +0.0822 ± 0.0006 for Jul2829, and (j)M = 
-0.0142 ± 0.0007 and (j)be = +0.073 ± 0.001 for 
JulSOB. BS enters eclipse later and reappears from 
eclipse earlier on JulSOB, indicating that the disk 
radius was smaller on that night. 

Each pair of {4>bi, (f'te) values maps into a single 
position in the orbital plane for an assumed binary 
geometry (i,g). The x-y positions correspond- 
ing to the measured {4>bi, (f'be) values are indicated 
by diamonds in the binary diagram shown in the 
right-hand panel of Fig. [S] The measured BS po- 
sitions consistently fall along the gas stream tra- 
jectory for the assumed mass ratio q — 0.19. The 



circles that pass through these points are Rbs — 
(0.480+0.015) Rli and Rbs = (0.S80+0.015) Rli, 
respectively for Jul2829 and JulSOB. The differ- 
ence between the two values is formally significant 
at the 6-CT confidence level. Thus, there is evidence 
that the disk shrank at outburst onset (i.e., the 
night before outburst maximum), in agreement 
with the expectations of MTIM. 

It shall be noticed that this result comes mainly 
from one light curve (JulSOB). Given the large 
flickering amplitude of V2051 Oph, one cannot ex- 
clude the possibility that the BS eclipse phases in 
this single light curve are affected by flickering. We 
attempted to strengthen the JulSOB result by re- 
peating the analysis for the V-band light curves of 
the same night. While the measured (j)bi is in good 
agreement with the B-band measurement, (j)be is 
unfortunately lost in the flickering on a slow egress 
slope. A comparison of the eclipse shape of the B- 
and V-band light curves of 2000 Jul SO indicates 
that the BS is less compact in the V-band than in 
the B-band. 

We further tested whether the reduction in 
brightness and radius of JulSOB was a common, 
short-lived effect caused by rapid changes in mass 
transfer rate (on timescales shorter than that of 
an outburst). We measured BS eclipse phases in 
the individual light curves of the previous nights, 
while the star was in quiescence. We took par- 
ticular attention to the first two eclipses on 2002 
Jul 28, for which the star had a brightness level 
comparable to that of the JulSOB light curve. In 
all cases the BS eclipse phases (and inferred disk 
radius) are consistent with the measurements ob- 
tained when we combine all data taken on 2000 Jul 
28-29. While the out-of-eclipse flux varies by ~ 20 
per cent among the individual light curves, we flnd 
no evidence of decrease in disk radius nor correla- 
tion of disk radius with out-of-eclipse brightness. 
This suggests that the observed reduction in disk 
radius reflects a change in disk structure particular 
of that night. 

4.2. Accretion disk structure 

4.2.1. The 2000 Jul Outburst 

B-band median light curves (dots with error 
bars) and eclipse mapping model curves (solid 
lines) are presented in the left panels of Fig. HI The 
right panels shows the corresponding eclipse maps 
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in a common logarithmic greyscale. The sequence 
of echpse maps ahows us to trace changes in the 
accretion disk structure at five different occasions 
along the outburst: in quiescence (Jul2829), and 
at outburst onset (JulSOB), maximum (JulSlB) 
and on two consecutive nights along the decline 
(AugOf and Aug02). 

The light curve and eclipse map in quiescence 
are very similar to those of Baptista & Borto- 
letto (2004), with enhanced emission along the gas 
stream trajectory inwards of the BS position. The 
tip of the asymmetric emission is consistent with 
the position where the ballistic stream trajectory 
intercepts the 0.48 Rli quiescent disk radius. The 
disk becomes fainter, smaller and largely asym- 
metric on Jul30B, with emission only along the 
gas stream trajectory. The lack of emission from 
the inner disk regions suggest a reduction of ac- 
cretion onto the WD, possibly as a consequence of 
redistribution of matter and angular momentum 
in the disk following outburst onset. 

The disk brightness distribution changes signif- 
icantly on the timescale of one day, from outburst 
onset (Jul30B) to outburst maximum (JulSlB). 
The wide, V-shaped and fairly symmetric eclipse 
of JulSlB maps into a broad brightness distri- 
bution that fills the primary Roche lobe, with 
a small asymmetry towards the trailing side of 
the disk (the one containing the stream trajec- 
tory). The brightness distribution is not concen- 
trated in the inner disk regions as it would be ex- 
pected for an opaque steady-state disk following 
the T{R) oc i?"3/4 (g g^^ 2002). 
The eclipse becomes progressively narrower and 
U-shaped along the following nights (AugOl and 
Aug02). The corresponding brightness distribu- 
tions show the inward cooling and fading of the 
outer disk regions whereas the inner disk remains 
at the same brightness level of outburst maximum. 
On Aug02 the outer disk is faint enough that 
an asymmetric and azimuthally extended emission 
from the BS region becomes clear. The evolution 
of the disk brightness distribution after outburst 
maximum is similar to that of the long-period 
dwarf nova EX Dra (Baptista & Catalan 2001). 

V-band light curves and eclipse maps of the 
2000 Jul outburst at outburst onset (JulSOV) and 
maximum (JulSlV) are shown in Fig. [5] The 
JulSlV light curve and eclipse map are similar 
to the corresponding B-band outburst maximum 
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Fig. 4. — B-band light curves and eclipse maps 
along the 2000 Jul outburst: in quiescence 
(Jul2829), at outburst onset (JulSOB), outburst 
maximum (JulSlB) and during decline (AugOl 
and Aug02). Left: median light curves (dots 
with error bars) and corresponding eclipse map- 
ping models (sohd lines). Vertical dotted lines 
mark the contact phases of the WD eclipse. Hori- 
zontal ticks depict the uneclipsed flux in each case. 
Right: the corresponding eclipse maps in a loga- 
rithmic greyscale common to all maps. Brighter 
regions are indicated in black, fainter regions in 
white. A cross marks the center of the disk; dot- 
ted lines show the primary Roche lobe, the gas 
stream trajectory for g = 0.19 and a reference disk 
of radius R\,s = 0.48 Rli] the secondary is to the 
right of each map. The horizontal bar indicates 
the logarithmic intensity level of the grey scale. A 
dashed contour line is overploted on each eclipse 
map to indicate the S-cr confidence level region. 



data. The JulSOV light curve has steeper ingress 
and egress slopes than the JulSOB light curve and 
the resulting eclipse map shows an additional light 
source at disk center aside of the enhanced emis- 
sion along the stream trajectory. We investigated 
whether this could be a consequence of underes- 
timating the contribution of the WD to the V- 
band light curve with the following exercise. We 
assumed progressively larger WD contributions to 
the JulSOV light curve and computed eclipse maps 
from the resulting WD-subtracted light curves (see 
Sect. 13. ip . The additional contribution at disk 
center is present in the eclipse map even at the 
limit where the WD-subtracted light curve starts 
to show a reverse slope at WD ingress/egress 
phases (signaling that too much flux was removed 
in the light curve decomposition process). Thus, 
we conclude that the extra light at disk center in 
the Jul30V light curve is not related to the WD. 
Since there are no V-band observations in the pre- 
vious nights, it is not possible to check whether the 
reduction in flux of the inner disk regions observed 
in the JulSOB B-band map is also present in the 
V-band. The tip of the emission pattern along the 
stream trajectory on JulSOV occurs at a radius 
smaller then the 0.48 Rli quiescent disk radius, 
in agreement with the inferred reduced radius at 
that date fSect. Oj) . 

4.2.2. The 2002 Aug outburst 

B-band light curves and corresponding eclipse 
maps along the 2002 Aug outburst are shown in 
Fig. m The greyscale of the eclipse maps is the 
same as in Fig. IH This outburst is poorly sampled 
by the amateur astronomers (Fig. [T]) and by our 
observations (Fig. [2]). It is possible to associate 
the Aug04 data with the quiescent, pre-outburst 
stage and the Aug07 data with the late decline 
stage of the short outburst. But it is not clear 
whether the Aug05 data corresponds to outburst 
maximum or not. Fig. [2] shows that there was 
no perceptible change in out-of-eclipse brightness 
over the > S hr long observing run of that night, in 
contrast with the clear decrease in brightness ob- 
served along the observations of the 2000 Jul out- 
burst. The constancy of the out-of-eclipse bright- 
ness suggests that the Aug05 observations framed 
a fairly stable brightness stage of the outburst pre- 
sumably not too far from outburst maximum. 

The symmetric and broad eclipse of AugG4 



maps into an axi-symmetric disk brightness distri- 
bution with no evidence of BS or enhanced emis- 
sion along the gas stream trajectory, in contrast 
with the 2000 Jul quiescent data. Also, the bright- 
ness distribution extends up to a smaller radius 
than in 2000 Jul. This eclipse map is reminiscent 
of that of the 'faint' quiescent state of Baptista & 
Bortoletto (2004), which may be attributed to a 
lower (long-term) mass transfer rate at that epoch. 
The Aug05 light curve shows an asymmetric V- 
shaped eclipse with two low-amplitude and ex- 
tended bulges at ingress and at egress. This leads 
to an eclipse map with an elongated annular struc- 
ture of asymmetric brightness distribution super- 
imposed on an axi-symmetric broad baseline dis- 
tribution. This disk brightness distribution bears 
some resemblance with the eclipse map of OY Car 
on the rise to maximum of a normal outburst (a 
ring-like structure, see Rutten et al. 1992b) and 
with the Hell A 4686 eclipse map of IP Peg at 
outburst maximum (a two-armed spiral structure 
on top of an extended brightness distribution, see 
Baptista et al. 2000). The disk radius increases 
with respect to Aug04 but is far from filling the 
primary Roche lobe. This seems in accordance 
with the markedly lower amplitude of this out- 
burst. The Aug07 light curve has a U-shaped sym- 
metric eclipse. The disk has decreased in size back 
to the pre-outburst radius but the inner disk re- 
gions remain at the intensity level of the Aug05 
map. The brightness distribution is skewed to- 
wards the LI point to account for an eclipse width 
at half-depth slightly wider than the eclipse of the 
white dwarf. 

4.3. Radial intensity distribution 

4.3.1. The 2000 Jul outburst 

A more quantitative description of the disk 
changes during outburst can be obtained by an- 
alyzing the evolution of the radial intensity dis- 
tribution. The left panels of Fig. [7| illustrate the 
evolution of the radial intensity distribution of the 
B-band maps along the outburst. We divided the 
eclipse maps in radial bins of 0.0Si?Li and com- 
puted the median intensity at each bin. These 
are shown as interconnected circles in Fig[71 The 
dashed lines show the ilcr limits on the average 
intensity. The labels are the same as in Fig. ID The 
large dispersion seen in the intermediate regions 
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Fig. 5. — V-band light curves and eclipse maps of 
the 2000 Jul outburst at outburst onset (Jul30V) 
and outburst maximum (JulSlV). The notation is 
similar to that of Fig. 21 
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Fig. 6. — B-band light curves and eclipse maps 
along the 2002 Aug outburst. The notation is sim- 
ilar and the greyscale is the same as that of Fig.[31 



of the JulSOB map and in the outer regions of the 
Aug02 map reflect the large asymmetries present 
in these eclipse maps (gas stream and bright spot, 
respectively). 

In order to quantify the changes in disk size 
during outburst we defined the outer disk radius 
in each map as the radial position at which the in- 
tensity distribution is logics — —4.8, which corre- 
sponds to the maximum intensity of the BS in the 
Jul2829 quiescence eclipse map. The computed 
outer disk radius is shown as a vertical tick mark 
in each panel. As a reference, the radial position 
of the quiescent BS is marked by a vertical dotted 
line. 

The disk shrinks from Rd = (0.48 ± 0.02) Rli 
in quiescence to Rd — (0.40 ± 0.11) at out- 
burst onset (JulSOB), underscoring the result of 
Sect. 14.11 The accretion disk then expands and 
reaches Rd = (0.78 ± 0.07) Rli at outburst max- 
imum (JulSlB). The disk decreases in radius to 
Rd = (0.70 ± 0.09) Rli on AugOl and remains at 
that radius on the following night (Aug02). 

The DIM predicts that heating and cooling 
wave-fronts propagate through the disk during 
the transitions between the low- viscosity quiescent 
state and the high-viscosity outburst state. The 
eclipse maps may be used to measure the move- 
ment of transitions fronts. In order to test for 
the presence and to trace the movement of tran- 
sitions fronts in the accretion disk, we defined an 
arbitrary reference intensity level log// = —4.77. 
This choice is justified by the following arguments. 
// > hs is required in order for it to trace the 
movement of regions inside the outer disk radius. 
Furthermore, If cannot be much higher than the 
chosen value, otherwise the whole of the bright- 
ness distribution on JulSOB would be below the 
reference intensity level. Choosing other // val- 
ues in the above range leads to results which are 
indistinguishable from the ones present here. 

In order to minimize the possible contribution 
of the BS and gas stream emission to the disk 
brightness distribution we calculated the symmet- 
ric disk-emission component. The symmetric com- 
ponent is obtained by slicing the disk into a set 
of radial bins and fitting a smooth spline func- 
tion to the median of the lower quartile of the 
intensities in each bin. The spline-fitted intensity 
in each annular section is taken as the symmet- 
ric disk-emission component. This procedure pre- 
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Fig. 7. — The evolution of the B-band radial in- 
tensity distribution through the 2000 Jul outburst. 
Left: Average intensity distributions. Dashed 
lines show the ilcr limits on the average intensity. 
A vertical dotted line indicates the radial position 
of the BS in quiescence, Rb, = OASRli (Ju12829), 
while vertical tick marks indicate the radius at 
which the distribution reaches the BS reference 
intensity level \og Its = —4.8. Right: the sym- 
metric disk-emission distribution, obtained from a 
cubic spline fit to the median of the lower quar- 
tile of the radial brightness distribution. Vertical 
tick marks indicate the radial position at which 
the disk intensity is log// ~ —4.77. 



serves the baseline of the radial distribution while 
removing all azimuthal structure. In doing this we 
are implicitly assuming that the global changes 
in the structure of the accretion disk with time 
(i.e., transition fronts and alike) are roughly axi- 
symmetric. The statistical uncertainties affecting 
the fitted intensities are estimated with the Monte 
Carlo procedure described in Sect. 13.21 

The radial position at which the intensity distri- 
bution falls below the reference intensity level //is 
indicated by vertical tick marks in the right panels 
of Fig. [71 It is seen that the radial position of the 
reference intensity changes significantly along the 
outburst, moving outwards on the rising branch 
and inwards along the decline. These changes sug- 
gest the propagation of an outward-moving heat- 
ing wave on the rise and of an inward-moving cool- 
ing wave on the decline. 

Assuming that the changes in the radial posi- 
tion of the reference intensity level represent the 
changes in position of the heating and the cooling 
waves, we used the measured positions together 
with the inferred time interval between consecu- 
tive eclipse maps to estimate the velocities of the 
waves along the outburst. The reference time as- 
sociated to each eclipse map is the average of the 
mid-eclipse times of all eclipses included in the cor- 
responding median light curve. The errors in the 
derived velocities have contributions from the un- 
certainties in the time interval between consecu- 
tive eclipse maps and the uncertainties in the mea- 
sured radial positions. The inferred speed of the 
heating wave is Vheat > (1.56 ± 0.05) kms~^. We 
quote this value as a lower limit because one can- 
not discard the possibility that outburst maximum 
(or largest disk radius) occurred earlier than our 
2000 Jul 31 observations. The inferred speed of the 
inward-moving cooling wave is Vcooii = (—0.30 ± 
0.07)fcms~i and Vcooi2 = (-0.90 ± 0.09) A:to s"!, 
respectively 1 and 2 nights after maximum. The 
results are summarized in Table [5] We observe an 
acceleration of the cooling wave as it travels across 
the disk, in contradiction with the predictions of 
the DIM (Menou et al. 1999). This result is sta- 
tistically significant at the 6-a confidence level. 

In terms of the a disk formulation of Shakura 
& Sunyaev (1973), the non-dimensional viscosity 
parameter of the high state, atot, can be written 
as the ratio of the speed at which the heating front 
travels across the disk, Vheat , and the sound speed 
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Fig. 8. — Left: radial intensity distributions 
through the 2002 outburst. The notation is the 
same as in Fig. [71 
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where Tj is the temperature of the heating front. 
Assuming Tf = 18000 iiT (Menou et al. 1999), we 
find ahot ^ 0.13, comparable to the measured vis- 
cosity parameter in quiescence acooi — 0.16 (Bap- 
tista & Bortoletto 2004). Unrealistically low Tf 
values (< 1000 K) are required in order to ob- 
tain ahot significantly larger than acooi- This re- 
sult is in contrast with the predictions of DIM, 
which requires that the viscosity in outburst be 
significantly larger than in quiescence (i.e., ahot ^ 
5 acooi)- 

4.3.2. The 2002 Aug outburst 

The radial intensity distributions along the 
2002 Aug outburst are shown in Fig. [51 We as- 
sumed the same reference intensity levels lbs and 
If as before to trace the changes in disk radius and 
to estimate the velocity of the heating and cooling 
waves. The derived Vf values (and their quoted 
errors) do not account for the uncertainty in as- 
sociating the Aug05 eclipse map with outburst 
maximum and, therefore, should be looked at with 



caution and skepticism. Moreover, there is a fur- 
ther uncertainty in the measured heating wave 
speed in this case because the Aug04 data does 
not corresponds to outburst onset but to a pre- 
outburst, quiescent state. The disk expands from 
Rbs = (0.35 ± 0.02) Rli in quiescence (Aug04) 
to Rbs = (0.50 ± 0.08) Rli at maximum observed 
brightness (Aug05) and shrinks back to the Aug04 
radius two nights after that (Aug07). We esti- 
mated an outward-moving heating wave speed of 
Vheat ~ 0.31 km and an inward-moving cooling 
wave speed of Vcooi ^ —0.15 km s^^. 

4.4. Radial temperature distribution 

A simple way to test theoretical disk models 
is to convert the intensities in the eclipse maps 
to blackbody brightness temperatures, which can 
then be compared to the radial run of the effec- 
tive temperature predicted by steady state, opti- 
cally thick disk models. It is important to bear 
in mind that B- and V-band brightness temper- 
atures are a good approximation to the effective 
temperature only for optically thick disk regions 
(see, e.g., Baptista et al. 1998b). Because it is 
hard to assess the optical depth of the disk gas 
along the outburst with the data at hand, the re- 
sults from this section should be looked upon with 
some caution. 

Figs. [Wroi show the evolution of the disk ra- 
dial temperature distribution along the 2000 Jul 
outburst in a logarithmic scale. The temperature 
distributions of the 2002 Aug outburst are shown 
in Fig. [m The blackbody brightness tempera- 
ture that reproduces the observed surface bright- 
ness at each pixel was calculated assuming dis- 
tances of 146 pc (Vrielmann et al. 2002a) and 
of 92 pc (Saito & Baptista 2006). These are re- 
spectively shown in the left- and right-hand panels 
of Figs. [Mm We neglected interstellar redden- 
ing, since there is no sign of interstellar absorp- 
tion feature at 2200A in spectra obtained with the 
Hubble Space Telescope (Baptista et al. 1998a). 
The disk was divided in radial bins of 0.03 
and a median brightness temperature was derived 
for each bin. These are shown as interconnected 
symbols in Figs. [MTTl The dashed lines show 
the ±1(T limits on the average temperatures. The 
larger tr- values in the Jul30B, Jul30V and Aug02 
distributions reflect the aziniuthal asymmetries 
caused by the gas stream and BS in these maps. 
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Fig. 9. — The evolution of the radial bright- 
ness temperature distribution along the 2000 
Jul outburst, assuming distances of 146 pc (left 
panel) and 92 pc (right) to the binary. Steady 
state disk models for mass accretion rates of 
M = 10-^ 10-10 and lO-^i Mej/r"! are plotted 
as dotted lines for comparison. A dot-dashed line 
marks the critical temperature above which the 
gas should remain in a steady, high mass accre- 
tion regime according to DIM. 



Steady-state disk models for mass accretion rates 
of M = 10-^10-l° and 10"" MQyr''^ are plot- 
ted as dotted lines for comparison. These models 
assume Mi = 0.78 M0 and a primary radius of 
Ri = 0.0103 i?0 (Baptista ct al. 1998a). 

The computed brightness temperatures depend 
on the assumed distance to the binary. For a 
given (fixed) observed flux, the disk is fainter (and 
cooler) if the distance is smaller. The tempera- 
tures and mass accretion rates quoted in the re- 
mainder of this section are for a 92 pc distance to 
the binary. 

Let's first discuss the 2000 Jul outburst dis- 
tributions. The quiescent disk (Jul2829) shows 
a flat temperature distribution for R < 0.3 Rli 
reminiscent of those seen in other dwarf novae 
(Wood et al. 1986,1989). The temperatures range 
from ~ 7000 K in the inner disk (i? = 0.1 Rli) to 
~ 5000 K in the outer disk regions {R = 0.4 Rli)- 
The whole disk brightens at outburst maximum 
(Jul31B and Jul31V), but the outer disk regions 
become relatively brighter and hotter, leading to a 
temperature distribution even flatter than in qui- 
escence for R < OA Rli- This is in marked con- 
trast with what is seen, e.g., in Z Cha (Horne 
& Cook 1985, Wood et al. 1986) and OY Car 
(Rutten et al. 1992b), the accretion disks of 
which transition from a flat brightness tempera- 
ture distribution in quiescence to a steep distri- 
bution in good agreement with the T oc R~^/'^ 
law in outburst. Along decline the V2051 Oph in- 
ner disk remains at about the same temperature 
(~ 10 000 K) while the flat portion of the distri- 
bution progressively recedes towards smaller radii 
as the outer disk regions cool down. Except for 
the inner disk regions {R < 0.2 Rli) - where the 
extra light source result in higher brightness tem- 
peratures - the V-band distributions are in agree- 
ment with their B-band counterparts, suggesting 
that, at least on 2000 Jul 30 and 31, the outer 
disk was optically thick and the derived brightness 
temperatures are a good approximation to the gas 
effective temperature. 

The temperature distributions share a com- 
mon feature; all show a flat inner portion which 
turns into a steep gradient closely following the 
T oc 7?""^/^ law in the outer disk regions. We fit- 
ted steady-state disk models to the steep regions 
of the temperature distributions to infer mass ac- 
cretion rates of M= (1.0 ± 0.1) x 10"^° Mgyr"^ 
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in quiescence, (6.4 ±0.1) x 10 Mq yr ^ at out- 
burst maximum, (3.7 ± 0.5) x 10"^" Mq yr~^ on 
AugOl, and (0.9±0.1) x 10"!° Mq yr'^ on Aug02. 
In the framework of the MTIM, the evolution of 
the outer disk temperature distribution could be 
interpreted as the response of a viscous disk to a 
sudden increase in mass transfer rate by a factor 
of 6.5. This corresponds to a brightness increase 
of AB — 2 mag, in line with the observed 2000 Jul 
outburst amplitude. 

Let's now turn our attention to the 2002 Aug 
outburst. The morphology of the temperature dis- 
tribution and its evolution with time arc similar to 
those of the 2000 Jul outburst, although the disk is 
cooler at maximum observed brightness (Aug05). 
The temperatures in the quiescent disk (Aug04) 
are comparable to those of Jul2829. The tem- 
perature distribution is flatter at maximum ob- 
served brightness, with T ~ 8 000K in the inner 
disk regions (flat portion of the distribution) and 
T ~ 5 700K at i? = 0.4 Rli- At the late de- 
cline stage (Aug07) the temperatures in the inner- 
most disk regions (i? < 0.15 Rli) are still com- 
parable to those at maximum brightness, while 
the outer disk has returned to the temperatures 
at quiescence. Fitting steady-state disk models 
to the steep outer disk regions of the tempera- 
ture distributions one finds mass accretion rates 
of M= (0.8 ± 0.1) X 10"^° Mq yr-^ in quiescence, 
(1.7±0.1) X 10"^" Mq yr~^ at maximum observed 
brightness, and (0.81 ± 0.05) x lO^i" Mq yr-^ at 
late decline. This could be interpreted as the re- 
sponse of a viscous disk to a sudden increase in 
mass transfer rate by a factor of 2.1, correspond- 
ing to a brightness increase of AB — 0.8 mag, 
again in line with the observed 2002 Aug outburst 
amplitude. 

Although the nominal temperatures and mass 
accretion rates change if one instead assumes the 
146 pc distance, the qualitative results remain the 
same. 

According to the DIM, there is a critical effec- 
tive temperature, Teff(crit), below which the disk 
gas should be while in quiescence in order to allow 
the thermal instability to set in, and above which 
the disk gas should remain in a steady, high mass 
accretion regime (e.g., Warner 1995). In other 
words, outbursting accretion disks should be hot- 
ter than Teff(crit). The Tefi(crit) relation is plot- 
ted as a dot-dashed line in each panel of Figs.lMTTI 
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Fig. 10.— The 2000 Jul radial brightness tem- 
perature distributions for the V-band data. The 
notation is the same as in Fig. [9l 
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Fig. 11. — Radial brightness temperature distri- 
butions for the 2002 Aug outburst. The notation 
is the same as in Fig. [51 
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The distance to the binary is the crucial pa- 
rameter for the Toff(crit) test, li d = 146 pc, the 
evolution of the radial temperature distribution is 
consistent with the DIM: the disk temperatures 
are everywhere below Tcff(crit) in quiescence and 
the outbursting disk shows temperatures clearly 
above that limit. However, for a 92 pc distance 
the outburst occurs at temperatures below (2002 
Aug) or barely at (2000 Jul) Tcff (crit). The Aug05 
outbursting disk temperatures are systematically 
lower than roff(crit) for distances d < 120 pc. 

4.5. The uneclipsed flux 

The uneclipsed flux is displayed as a horizon- 
tal tick centered at phase zero in the left panels of 
Figs. mis] and m The measured values are listed in 
Table m We estimated the fractional contribution 
of the uneclipsed component to the total flux by 
computing the ratio of the uneclipsed flux, F^n, to 
the average out-of-eclipse flux level of the corre- 
sponding light curve, Fout- The results are shown 
in the third column of Table [H 

The uneclipsed component accounts for a small 
fraction (1.5 percent) of the total system bright- 
ness on Jul2829 but increases to 13 percent on 
JulSOB, the highest relative contribution mea- 
sured. This is in line with the significant changes 
in disk radius and structure observed on JulSOB 
and corroborates the idea that this night corre- 
sponds to outburst onset. Although the abso- 
lute value of Fun is larger at outburst maximum, 
its fractional contribution decreases steadily since 
2000 Jul 30 and becomes neghgible on 2000 Aug 
2. The behaviour in the V-band is analogous 
to that in the B-band, with a non-negligible Fun 
value on Jul30V increasing to a maximum (in ab- 
solute terms) on the following night. The un- 
eclipsed component is negligible throughout the 
fainter 2002 Aug outburst. We draw attention to 
the fact that the increase in Fun on the 2000 data 
precedes outburst maximum by about one day. 

The fact that the uneclipsed component is neg- 
ligible in quiescence and its variability during out- 
burst indicate that this light does not arise from 
the secondary star. Significant (and variable) 
uneclipsed components during outburst were re- 
ported for OY Car (Rutten et al. 1992b) and 
EX Dra (Baptista & Catalan 2001). Two pos- 
sible explanations to account for the observed be- 
haviour would be a variable and largely uneclipsed 



disk wind or a flaring of the accretion disk during 
outburst. 

Detailed simulations by Wood (1994) show that 
eclipse maps obtained including an uneclipsed 
component or assuming a flared disk may lead 
to equally good fits to the data light curve, mak- 
ing it hard to distinguish between the two models. 
Further simulations by Baptista & Catalan (2001) 
show that modeling the eclipse of a flared disk 
with the assumption of a flat disk leads to the ap- 
pearance of a spurious uneclipsed component, the 
flux of which scales with the disk opening angle. 

For a mass accretion rate of ~ lO~^°M0j/r~^ 
(sect. I4.4p . we estimate a quiescent disk opening 
angle of 0d =i 2° (Meyer & Meyer-Hofmeister 1982; 
Smak 1992). To be accounted for by a flaring of 
the accretion disk, the observed increase in un- 
eclipsed flux by a factor ~ 7 from Jul2829 to 
Jul30B would demand an increase in 9d by fac- 
tors of 7-25 (see Baptista & Catalan 2001), leading 
to obscuration of a large portion of the accretion 
disk, in particular the disk center. This is hard 
to reconcile with the fact that the white dwarf at 
disk center was still visible in the 2000 Jul 30 light 
curves. Thus, an outburst-driven disk wind seems 
the most plausible explanation for the observed 
uneclipsed flux. 

5. Discussion 

5.1. The distance to the binary 

Because the radial temperature distribution 
and the consequent comparison of disk temper- 
atures with Tcff (crit) depend on the assumed dis- 
tance to the binary, we now turn our attention to 
the distance estimates in the literature. 

Watts et al. (1986) found a distance between 
90 and 150 pc assuming a steady-state optically 
thick model for the quiescent accretion disk of 
V2051 Oph. However, the UV-optical spectrum 
of V2051 Oph is dominated by strong emission 
lines and a Balmer jump in emission (Watts et al. 
1986, Baptista et al. 1998a) ~ indicating impor- 
tant contribution from optically thin emitting re- 
gions ~ and the flat temperature distribution of 
the inner disk regions in quiescence is in clear dis- 
agreement with the T oc i?"^/'' law of steady-state 
opaque disks (e.g.. Fig. [9|). Berriman et al. (1986) 
compared the eclipse shape in the H-band with 
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computed models to find distances in the range 
120-170 pc. They assumed that the red dwarf 
has a flux density of 1.5 mJy at the H-band and 
that cool stars have approximately constant sur- 
face brightness at the H-band. Their models are 
a poor match to the observed eclipse shape. Both 
distance ranges are affected by arguable assump- 
tions and/or too simple models. 

Vriclmann et al. (2002a) applied a physical pa- 
rameter eclipse mapping method (PPEM) to UB- 
VRI light curves of V2051 Oph in order to derive 
the spatial distribution of temperature and surface 
density in its accretion disk and to find a distance 
of d = 146 ± 20 pc to the binary. This method re- 
quires the assumption of an a priori spectral model 
for the disk emission relating the temperature and 
surface density to the observed surface brightness 
in the UBVRI passbands. The adopted spectral 
model is an isothermal, pure hydrogen slab of gas 
in LTE including only bound-free and free-free H 
and H~ emission. This model is not adequate to 
handle emission from disk atmospheres with verti- 
cal temperature gradients. In particular, it might 
fail to model a temperature inversion in the disk 
photosphere (i.e., a hot disk chromosphere with 
strong emission lines), as it seems to be the case 
in V2051 Oph (Vrielmann et al. 2002a; Saito & 
Baptista 2006) and other quiescent dwarf novae. 
The application of this method and model to the 
dwarf nova HT Gas leads to a distance 50 per 
cent larger then that inferred from the Nal and 
TiO absorption spectrum of the secondary star 
(Vrielmann et al. 2002b). The PPEM distance to 
V2051 Oph may be overestimated as well. More- 
over, Vriclmann et al. (2002a) mistakenly adopted 
a WD radius of i?i = 0.0244 Rq, a factor of 2.4 
larger than the i?i = 0.0103 Rq value quoted by 
Baptista et al. (1998a). As a result, the WD 
ingress and egress features in the light curve be- 
come artificially longer and the PPEM code asso- 
ciates more surface flux to the WD than it should. 
The WD accounts for a significant fraction of the 
V2051 Oph optical light in quiescence (see Fig. 3 
of Vrielmann et al. 2002a) and dominates the pro- 
cedure used to estimate the distance. The extra 
WD surface flux has to be compensated for by in- 
creasing the distance to the binary. 

We are therefore left with the d — 9213° pc dis- 
tance estimate of Saito & Baptista (2006). This 
was obtained by fitting white dwarf atmosphere 



models of solar composition to the extracted UV- 
optical WD spectra at the time the WD was most 
clearly seen in V2051 Oph (Baptista et al. 1998a; 
Saito & Baptista 2006). In the absence of a more 
robust determination we take this as the best 
available distance estimate to V2051 Oph. 

5.2. The cause of the V2051 Oph out- 
bursts 

Eclipse mapping analysis along two outburst 
cycles in V2051 Oph allows us to test the pre- 
dictions of dwarf novae outburst models listed in 
Sect.[T] Here we discuss the cause of the outbursts 
of V2051 Oph under the light of the results pre- 
sented in the previous sections. 

MTIM predicts that the disk shrinks at out- 
burst onset (e.g., Warner 1995 and references 
therein). We found that the disk decreased in 
brightness and radius at the night before the 2000 
Jul outburst maximum while the uneclipsed light 
in the B-band increased (from a negligible contri- 
bution) to 13 per cent of the total light. It is hard 
to avoid the interpretation that the changes in disk 
structure and uneclipsed flux are a consequence of 
the outburst onset and the following conclusion 
that our result confirms the MTIM prediction. 

DIM requires that auot ^ 5 ol^ooI (e.g., Lasota 
2001). In MTIM, there is no reason for the viscos- 
ity parameter to change from quiescence to out- 
burst. We measure othot — olcooI — 0.13, confirm- 
ing another MTIM prediction. In favor of DIM, 
one could argue that the true Vheat may be much 
larger than measured, contributing to increase 
ahot relative to acooi- However, there is no ob- 
servational support for an ahot significantly larger 
than inferred. Vrielmann & Offutt (2003) mea- 
sured a heating wave speed of Vheat = +1.8 km s^^ 
during a superoutburst in V2051 Oph, in good 
agreement with our result. One could alterna- 
tively dismiss the acooi measurement of Baptista & 
Bortoletto (2004) , leaving room for a much smaller 
<Xcooi value. Nevertheless, the good agreement be- 
tween the radial temperature distribution of the 
outer disk and the T oc R^^^^ law suggests that 
the outer regions of the V2051 Oph accretion disk 
are in a steady-state both in outburst and in qui- 
escence, implying that the disk viscosity there is 
always high. 

DIM predicts that the cooling wave decelerates 
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as it travels towards disk center (Menou et al. 
1999). In MTIM, the measured Vcooi reflects the 
local viscous timescale and may increase or de- 
crease depending on the radial dependence of the 
viscosity parameter a. We find that the cooling 
wave accelerates as it travels towards disk center, 
in clear disagreement with the DIM prediction. 

DIM predicts that an outbursting accretion 
disk is hotter than Teff(crit) (e.g., Warner 1995 
and references therein). There is no quies- 
cence/outburst temperature restriction in MTIM. 
For an assumed distance of 92 pc, we find that the 
whole of the low-amplitude 2002 Aug outburst oc- 
curs at temperatures below Teff(crit), and that the 
hottest parts of the accretion disk during the nor- 
mal 2000 Jul outburst barely reach Teff(crit). The 
2002 Aug outburst fails to meet the T > reff(crit) 
requirement for distances d < 120 pc. The mor- 
phology and time evolution of the radial temper- 
ature is similar on both outbursts, and there is 
no reason to believe they are powered by different 
physical mechanisms. If one admits the 2002 Aug 
outburst is not driven by a thermal-viscous disk 
instability, the same conclusion must be applied to 
the normal 2000 Jul outburst. The inconsistency 
between DIM and the observations here could be 
alleviated if one argues that the V2051 Oph ac- 
cretion disk is not optically thick during outbursts 
and that the derived brightness temperatures un- 
derestimate the disk gas effective temperatures. 

If taken separately, none of the individual re- 
sults would be enough to discard DIM as a pow- 
ering mechanism for the outbursts of V2051 Oph. 
Nevertheless, when taken in combination the 
above results make a strong case against DIM and 
in favour of MTIM as the cause of the outbursts 
in V2051 Oph. All evidences discussed here lead 
to the conclusion that the outbursts of V2051 Oph 
are driven by episodes of enhanced mass-transfer 
from its secondary star. 

Before one claims that V2051 Oph cannot be 
considered a 'true' dwarf nova, we remind it is an 
SU UMa star, showing not only 2 mag amplitude 
normal outbursts, but also longer and brighter su- 
peroutbursts (Kiyota & Kato 1998; Vrielmann & 
Offutt 2003). It is also worth noting that finding a 
dwarf nova the outbursts of which arc not powered 
by disk instabilities does not necessarily kills DIM, 
but warns that the current view about dwarf nova 
outbursts must be revised to include the possibil- 



ity that both outburst mechanisms coexist, per- 
haps on different dwarf nova subtypes that have 
not been clearly distinguished up to this point. 

We propose that the viscosity in the accretion 
disk of V2051 Oph is always high, independent of 
the disk mass inflow rate. In this scenario, there 
is no room for matter to accumulate in the outer 
disk regions nor for thermal-viscous instabilities 
to set in. Because of the low-mass transfer rate 
and the efficient accretion engine, the disk has low 
densities. A low density quiescent disk enables 
the (denser) infalling gas to 'overflow' the outer 
disk regions, creating a bright gas stream along 
the ballistic trajectory ahead of the disk rim (cf., 
Baptista & Bortoletto 2004) and making the trail- 
ing Imic of the disk significantly brighter than the 
leading lune - as observed by Watts et al. (1986). 
This is supported by numerical simulations of ac- 
cretion discs, which show that when the infalling 
gas stream is significantly denser than the initial 
disc gas, the stream penetrates the disc - allowing 
matter to be deposited at the inner disc regions 
- and no bright spot forms at disk rim (Bisikalo 
et al. 1998a,b; Makita et al. 2000). Because the 
energy released at the disk rim is small in compar- 
ison to that released along the stream closer to the 
WD, there is no pronounced orbital hump caused 
by anisotropic BS emission in the light curve - in 
contrast with the conspicuous orbital hump seen 
in the light curves of the low- viscosity disk dwarf 
novae OY Car (e.g.. Wood et al. 1989), Z Cha 
(e.g.. Wood et al. 1986) or IP Peg (e.g.. Wood 
& Crawford 1986). Furthermore, without strong 
anisotropic BS emission, there will be no observ- 
able increase in BS luminosity to signal a sudden 
enhancement of mass transfer rate at outburst on- 
set. Since the transferred matter is allowed to 
flow all the way down to the inner disk regions 
(i.e., at the circularization radius), a burst of en- 
hanced mass-transfer from the secondary star can 
lead to inside-out outbursts similar to that of 2000 
Jul. Variability in the mass transfer from the sec- 
ondary star is therefore responsible not only for 
the strong flickering and the long-term brightness 
changes (Baptista & Bortoletto 2004) but also for 
the infrequent normal (increase in M2 by factors 
5-10) and low-amplitude (increase in M2 by fac- 
tors 2-3) outbursts. Going one step further, there 
is no reason to believe the disk would switch to 
a low-viscosity state just before a superout burst. 
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Therefore, this scenario impHes that the rare su- 
perout bursts of V2051 Oph are also powered by 
(stronger) bursts of enhanced mass transfer from 
its secondary star. 

The increase in the luminosity of the BS at out- 
burst onset and the outburst type issue (inside-out 
versus outside-in) are no valid tests to distinguish 
between DIM and MTIM in dwarf novae such as 
V2051 Oph. Are there other, useful tests to distin- 
guish between the two outburst models? Section 1 
gives a list of observational tests for which the pre- 
dictions of both models are quite different (and, 
therefore, relatively c^asy to distinguish). Unfortu- 
nately, they are only applicable to the restricted 
group of eclipsing dwarf novae. Comparing the 
disc radius in quiescence and at outburst onset 
would be the most straightforward test, but also 
the hardest to perform because it requires observa- 
tions framing few hours around the highly unpre- 
dictable start of the outburst. The comparison of 
the disc brightness temperature distribution with 
the roff(crit) relation requires good knowledge of 
the distance to the binary and depends on the 
assumption that the brightness temperature is a 
good approximation to the disc gas effective tem- 
perature (or, alternatively, on the availability of 
an adequate model relating both temperatures). 
Tracing changes in speed of the cooling wave along 
decline is probably the less demanding test from 
the observational point of view - one needs good 
sampling of the decline branch of the outburst 
- but may be inconclusive if one finds that the 
cooling wave decelerates. Because ahot is easily 
inferred from the timescale of the outburst, the 
key aspect of the comparison of ahot and acool 
is the measurement of the viscosity parameter in 
quiescence. A promising step in this regard may 
be to perform spatially-resolved flickering studies 
of quiescent dwarf novae. The identification of a 
disc flickering component allows an estimate of the 
magnitude and the radial dependency of the disc 
viscosity (e.g., Baptista & Bortoletto 2004). 

The present results raise a set of exciting and 
challenging questions: How many other high- 
viscosity disk dwarf novae similar to V2G51 Oph 
are there? What fraction of the dwarf novae class 
do they account for? Why do their accretion 
disks have permanent high viscosity? Why and 
how the secondary stars of CVs change their mass 
transfer rates by large amounts on short, 1-2 days 



timescales? 
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Table 1 
Journal of the observations 



date 


HJD start 
(2.450.000+) 


At 

(b) 


band 


E ^ 
(cycle) 


phase 
range 




Quality ^ 


brightness 
state 


2000 Jul 28 


1754. 44d4 


5 


TO 

a 


136293 


A 1 A 1 A r A 

—0.19, +0.50 


726 


A 


quiescence 




1754.4894 


5 


TO 

a 


1 o/^on A 

ioDzy4 


A r A 1 A rr A 

—0.50, +0.50 


A /I 1 

941 


A 

A 


quiescence 




1754.5518 


5 


TO 

a 


135295 


A A 1 A C A 

— U.5U, +U.5U 


1 AT /I 

1074 


A 

A 


quiescence 









L) 




— U.ou, -rU.Zo 


Ooo 


A 


quiescence 


2000 Jul 29 


1755.3952 


5 


B 


1 o o o 

136308 


1 A A 1 (ATA 

+0.01, +0.50 


281 


A 

A 


quiescence 




1755.4257 


5 


B 


"1 o o r\ A 

136309 


ATA 1 A P" A 

—0.50, +0.50 


1079 


A 

A 


quiescence 




i ( 00.4OOZ 





TO 
L) 


ioDoiU 


— u.ou, +U.OU 


iU ( o 




quiescence 






c; 
o 


T{ 


1 ^fim 1 

-LOUO ± X 




8fin 

ouu 


A 


quicoccliuc 




1 755 fil 7fi 


K 
o 




1 3631 9 


— 43 +0 36 


848 


A 


mm PQPPTI pp 


9nnn Tni "^n 


L 1 OK) .'±±\JU 


o 


_L) 






you 


R 


/"(inpcr'PTi/^p 

quicoccliuc 




1756.5029 


5 


V 


136326 


-0.24, +0.50 


793 


B 


quiescence 




1756.5494 


5 


V 


136327 


-0.50, +0.30 


836 


B 


quiescence 




1756.6285 


5 


V 


136328 


-0.23, +0.45 


850 


C 


quiescence 


2000 Jul 31 


1757.4083 


5 


V 


136341 


-0.74, +0.50 


1296 


C 


outburst 




1757.4859 


5 


V 


(136342) 


-0.50,-0.09 


376 


B 


outburst 




1757.5116 


5 


B 


136342 


-0.09, +0.50 


627 


A 


outburst 




1 / / .04oo 





13 

_t> 


ioDo4o 


n _i_n 
— U.ou, H-U.O ( 


1 1 Qc: 
iioO 


A 


outburst 




1757.6152 


5 


V 


136344 


-0.43, +0.45 


949 


A 


outburst 


2000 Aug 1 


1758.4741 


5 


B 


136358 


-0.67, +0.50 


420 


B 


decline 




1758.5472 


10 


B 


136359 


-0.49, +0.50 


254 


C 


decline 




1758.6096 


15 


B 


136360 


-0.50, +0.14 


139 


C 


decline 


2000 Aug 2 


1759.4485 


5 


B 


136373 


-0.06, +0.50 


202 


A 


decline 




1759.4837 


5 


B 


136374 


-0.50, +0.50 


323 


B 


decline 




1759.5461 


5 


B 


136375 


-0.50, +0.10 


217 


A 


decline 


2002 Aug 4 


2491.4401 


10 


B 


148099 


-0.66, +0.50 


626 


B 


quiescence 




2491.5126 


10 


B 


148100 


-0.50, +0.75 


672 


B 


quiescence 


2002 Aug 5 


2492.4126 


10 


B 


148114 


-0.08, +0.50 


315 


A 


outburst 




2492.4491 


10 


B 


(148115) 


-0.50, +0.25 


255 


C 


outburst 




2492.5278 


10 


B 


148116 


-0.24, +0.16 


174 


C 


outburst 


2002 Aug 7 


2494.4529 


20 


B 


148147 


-0.40, +0.58 


266 


A 


decline 



^with respect to the ephemeris of eq. ([T]) . 

^sky conditions: A= photometric (main comparison stable), B= good (some sky variations), C= poor 
(large variations and/or clouds). 
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Table 2 

Measuring the speed of transitions waves 



eclipse map 


Rf/R 


Ll 


ARf/RLi 


At 


Vf 


sequence 








(days) 


{kms ^) 


Ju130Bh 


■^Jul 31B 


0.22 ^ 


0.72 


+0.50 


1.09 


+1.56 + 0.05 


Ju131Bh 


-> Aug 01 


0.72 ^ 


0.63 


-0.09 


1.03 


-0.30 + 0.07 


AugOl H 


Aug 02 


0.63 ^^ 


0.38 


-0.25 


0.94 


-0.90 + 0.09 


Aug04 h 


Aug 05 


0.34 ^ 


0.43 


+0.09 


0.97 


+0.31 + 0.04 


Aug05 H 


-> Aug 07 


0.43 


0.34 


-0.09 


2.00 


-0.15 + 0.02 



Note. — Rf is the radius at which the intensity falls below log// = 
—4.77; At is the time interval between the two consecutive eclipse maps; 
Vf = {ARf/At) is the speed of the moving wave (in kms~^). 



Table 3 
The uneclipsed flux 



Map J^un (mJy) Fun/-Fout 



Jul2829 


0.03 + 0.02 


0.015 


Jul30B 


0.20 + 0.04 


0.13 


Jul31B 


0.42 + 0.07 


0.04 


AugOl 


0.17 + 0.06 


0.02 


Aug02 


0.00 + 0.01 




Jul30V 


0.10 + 0.02 


0.04 


Jul31V 


0.62 + 0.06 


0.04 


Aug04 


0.01 + 0.01 




Aug05 


0.02 + 0.01 




Aug07 


0.00 + 0.01 
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